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| Abstract 



The parquet approximation in the matrix Higgs model is consid- 
ered. We demonstrate analytically that in the large N limit the par- 
rel 1 quet approximation gives an satisfying agreement with the exact re- 
sults. It is shown that the parquet planar series can be derived by 
means of the generating functional. 



*h : 1 Introduction 



The parquet approximation was proposed by Landau, Abrikosov and Kha- 
latnikov in order to develop a self-consistent method of studying the non- 
perturbative domain in quantum electrodynamics [I]. Later on, this ap- 
proximation has been used for various models of quantum field theory [2]. 
The parquet approximation leads to a closed system of integro-differential 
equations which have meaning not only for small but also for large values 
of coupling constant. The main problem that prevents direct application of 
the parquet approximation to an arbitrary gauge theory is that the parquet 
approximation violates the gauge invariance of the theory. 

* iskOziyatdin . phys . msu . su 
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Another non-perturbative approach used in the quantum field theory is 
the planar approximation, a method of studying theories like SU(N) QCD 
in the large N limit [3]. It enables us to understand most specific features 
of QCD jS], and unlike the parquet approach it does not break the gauge 
invariance. But to perform analytical investigations one has to calculate the 
sum of all planar diagrams, and this problem is solved only for a few simple 
models such as zero- and one-dimensional matrix model and two-dimensional 
QCD 0, H. 

There is a hope to construct a method that will combine advantages of 
both planar and parquet approaches. 

In jU] the planar parquet approximation was defined and applied to zero- 
dimensional matrix models with cubic and quartic interaction terms. It 
was demonstrated that the planar parquet approximation gives an excel- 
lent agreement with the exact results. The aim of this paper is to find out 
whether this approximation can give any sensible results in the matrix Higgs 
model. In order to answer this question we calculate Green functions within 
both approaches and compare them. 

The paper is organized as follows. In Section 2 we define the planar par- 
quet approximation and use it for studying the matrix Higgs model. In Sec- 
tion 3 the same model is considered by means of the steepest descent method. 
In Section 4 the generating functional for the whole parquet planar series is 
constructed. It is shown that this functional can be regarded as a restricted 
planar one. 



2 Planar parquet approximation 

In this section we define and investigate the planar parquet approximation 
for the zero- dimensional matrix Higgs model. To introduce it carefully it is 
useful to recall what is already known. 

The parquet planar approximation for the (i-dimensional matrix model 
with cubic and quartic interaction has been defined in j^J. Let us consider 
the matrix model with the action 

S = f d d x Tr ( -DM ■ dM + -m 2 M 2 ^=M 3 + — M 4 ) , (2.1) 

J \2 2 3^ 4iV J' K J 

here M is a hermitian matrix N x N. 
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The planar Green functions 



1 jDM^{M{ Xl )...M{x n ))eMS) 
ll n {x u ...,x n ) - mn ]Vl+n/2 /DMexp(-5) 1 j 

satisfy the planar Schwinger-Dyson equations [TUl HI] 

(- A +m 2 ) xi n n (xi, ... 

- An n+ i(si, xi,x 2 , ...,x n ) + gU- n+2 {x 1: xi, xi, x 2 , x n ) = 

n 

= - x m )n m _ 2 (x 2 , ...,x m _i)Il n _ m (a; m+ i, ...,x n ) = 0. (2.3) 

m=2 

This set of equations is infinite. So the purpose of the planar parquet ap- 
proximation is to find conditions when (J2.3|) becomes closed, i.e. the planar 
parquet approximation is defined as an approximative perturbative solution 
of ()2.3|) that takes into account only a part of the full series of coupling 
constant. 

It can be done by using a so-called skeleton expansion which contains only 
a subset of all planar diagrams. More precisely, the diagrams from this subset 
contains no bare propagators, three- and four-vertices insertions. The basic 
Green functions, i.e. the two-, three- and four-point functions are defined 
as solutions of a set of integro-differential equations. The zero-dimensional 
case reduces it to the set of algebraic equations. One can represent them 
graphically 



+ 



+ 





+ 
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(all contributions of tadpole diagrams are dropped out). Here the thick and 
thin lines represent the full (within the planar parquet approximation) and 
bare propagators respectively 

= D, = 1, 



are the full three- and four-point vertices, 



I 



// • I = V 



are the parts of the four-point vertex function that are 2PR in the t-channel 
(s-channel) and not 2PR in the s-channel (t-channel). The vertices V and 
H are related by the cyclic permutation of external points. 

Let us apply the foregoing considerations to the zero-dimensional matrix 
Higgs model with the action 

S = Tr (~\m 2 + j^M 4 ] . (2.4) 



The classical vacua are 

5S 



M = 0, 

- t N 2 
5M 1 Tr M 2 = -—. 

4g 
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There exist two possible ways to get the perturbative solution: we can 
write down planar parquet equations in the vicinity of the false vacuum 
(M = 0) and near the true vacuum ( Tr M 2 = N 2 /4g). 

In the first case we have the following set of equations on D, T^, H, V 



D = -1 -8gD 2 -AgD 4 T 4 
T 4 = -4g + H + V 
H = -AgD 2 T 4 + VD 2 T A 
V = -AgD 2 Y 4 + HD 2 T 4 



(2.5) 



As one can see this is a set of four equations for four variables. Excluding 
T 4 , V, H one get the following equation 



Ug A D b + \Qg 2 D b + \Ylg 2 D 4 + 20gD d + (1 + 20g)D 2 + 2D + 1 = 0. (2.6) 

It can be solved in the limit of small g. There exist two roots that have 
no singularities as g — ► 



£>(2) 



l2g + o(g), 
8g + o(g). 



Moreover, there exists a solution behaving like — as g 

9 

write down the following equation for a 



64a b + 16a 5 = 0, 



(2.7) 
0. One can 

(2.8) 



wherefrom a 



-- and 
4 



(2.9) 



The following table contains results of numerical calculations of ()2.6|) . 
The solutions exist if g < 0, 038. 



9 





i(r y 


10" 5 


10~ 4 


10~ 3 


10~ 2 


io- 1 


£>(D 


-1 


-1,000012 


-1,000120 


-1,001203 


-1,012331 


-1,169973 




D (2) 


-1 


-1,000008 


-1,000080 


-1,000801 


-1,008113 


-1,093576 




£>(3) 




-249997,0 


-24996,9 


-2496,9 


-246,9 


-21,8 
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The second case is more delicate. Consider a shift to the true vacuum 



M, 



a(3 




Ra/3 + \ -T-Iaf3, 



Tr R = 0. 



Hence, we get the following action 
S = I i R 2 



N 



For convenience one can rescale the fields R 



s - Tr ( — 4 

2 v^v^V 



N~~ 16g 
Q_ 

i.v"" T% 



(2.10) 



(2.11) 



It contains both cubic and quartic interaction terms. 

Thus, in the zero-dimensional case the planar parquet equations on D, 
T 3 , rV look like 



z? = i — — ^— ^ 3 r 3 - 2^ 2 - gD*T 4 



r, = 



2 

3^ 



r 4 = + if + V + DT; 



v ^ +jD 2 Vr3+jD 3 r 3_^ r3jD 2 
~4r4 



C2-12) 



y V = D 2 f 4 H + D 3 f \H + D 3 f 4 f § - #f) 2 f 4 - #f) 3 f § 
This set can be solved as g — > 

5 = l + ^ + o((?). 

This set of equations gives a sensible physical solution if g < 0, 037. 



(2.13) 
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1(T B 


10~ 5 


io~ 4 


io- a 


10~ 2 


io- 1 


D 


1 


1,000002 


1,000060 


1,000250 


1,002504 


1,030238 


-2,220680 



Hence, the true value of the propagator is 

D = < TrM 2 >=^ + i < Tr Q 2 >=^ + l + ^g + o(g). (2.14) 
Ag 2 Ag 2 4 



iHere L> = < Tr Q 2 >, T 3 = < Tr Q 3 >, T 4 = < Tr Q 4 > are Green functions full 
within the planar parquet approximation. 
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3 Planar approximation 

In this section we study the zero-dimensional matrix Higgs model by means 
of the approach proposed in A specific feature of this model is that it 
contains two-cut solutions parameterized by a special parameter |Z1 IH1 E] • 
In [4] one considered the large N limit in the model 

S= Tr Qm 2 + -|m 4 ^ . (3.15) 

By means of the steepest descent method the vacuum energy 

e -N*E (g) = Um f DMe S / 31 g) 
N-too J 

and Green functions 

1 JPMTr(M( Xl )...M M )eMS) 
G n { Xl ,...,x n ) - ^im /DMexp(-5) ld ?J 

were calculated. 

The results are presented below 

E (g) = J dXu(X) ^X 2 + gX 4 ^j- J J dfid\u(X)u(fi)\n\\ - (3.18) 

G 2p = [ X 2p u(X)d\, (3.19) 
here one introduced the eigenvalue density function u(X) such that 

u(X)dX = 1, 



x 

where X is the support of u(X). This density function is the solution of the 
following singular equation 

iA + 2sA 3 = v.p. [ ^-d/*. 
2 JxX-^i 

The model ()3.15|) admits so called one-cut solutions, i.e. when X has the 
form of a single segment (2a, 26). This support is uniquely defined by a set 
of algebraic equations. 



7 



The matrix Higgs model 

S = Tr (-\m 2 + j^M 4 ^ . (3.20) 

is interesting because it is the simplest model where exist both one-cut and 
multi-cut (two-cut) solutions. As it was mentioned in [SJ E] in the case of 
two-cut solution there exists certain freedom, to fix it one has to introduce 
an extra parameter. It is associated with the order parameter which governs 
the phase structure of the system. 

Consider the one-cut solution, i.e. let X = (2a, 26) = (t — s, t + s). There 
exist two solutions. The first one corresponds to symmetric support t — 0. 
It gives the density function 

u(X) = ±- (8a 2 g + AgX 2 - l) ^/{2af - A 2 , g > 0, (3.21) 
where a is subject to 

12^a 4 - a 2 - 1 = 0, (3.22) 
and the two-point Green function 

f 2b a 2 (4 + a 2 ) 1 1 

G 2 = dX X 2 u(X) = 1 ; a> = 7^ + ^ + 1-8^ + o(g). (3.23) 
J-2a 3 432# 2 Qg 

The second solution corresponds to the non-symmetric support 

WOO? 0<9< 1/60 (3 24) 

20g 15$ ' y - i \ ) 

The density function looks like 

u(X) = — (4gX 2 + AgtX + Agt 2 + 2gs 2 - l) V2tA - A 2 + s 2 - t 2 (3.25) 
and the two-point Green function is 

G 2 = - ^ + 3g*V + 3gsH 2 + ^ = 1 - 1 - Wg + o{g). (3.26) 

4 Id 4 4^ 

As in the previous section one can perform these calculations in the dif- 
ferent way: one can make the shift (j2.10|) to the true vacuum and apply the 
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same technique to the action (J2.ll)) . But this case can be solved only in the 
limit of small g. The results are presented below: 
the support X = (2a, 26) is 

/ 3 15 ( 3 - 27 ) 

the density function is 

u(\) = ^(2g\ 2 + 2g(a+b)\ + 3^g\ + 2g(a+b) 2 + 3 y /g(a+b)+g( y a-b) 2 + l)x 

x y/(2a-X)(X-2b), (3.28) 



7T 



and the behavior of the two-point Green function as g — > 0, taking into 
account 1)2. 14J) . is 

G2 = ^ + \ + Sg + o(g). (3.29) 

As one can observe two different values for G 2 (13.26)) and ()3.29|) are not 
the same, they show similar behavior as g — > 0. 

Consider the two-cut solution. The simplest case is the symmetric support 
X = (—26, —2a) U (2a, 26). By means of the previous procedure one can get 
the following: 
the support X is 

a = — £— , 6 = — =— , < 5- < 1/16, (3.30) 

the density function is 

u(A) = j^v 7 ^ 2 - !6^ 2 A 4 - 1 + I60, (3.31) 

27T 

and the two-point Green function 

G 2 = 1 (3.32) 

One can see that in this case the final result is exact. 

One has to notice that ()3~26jl . ()3~29j) and ()3~32jl have the same behavior 
as g — > 0. 
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4 Parquet planar generating functional 

1. Planar generating functional 

For the sake of simplicity we will study the model with the following 
action 

S = -Tr M 2 + -|-Tr M 4 . (4.33) 
2 47V v ' 

The planar Schwinger-Dyson equations in this case are 

n-2 

n n + g u n+2 = nin n _,_ 2 , n > 2. (4.34) 

i=0 

To solve them one introduces the following functional 

00 

F(x) = ^yn n . 

n=0 

where n = 1 and n 2fc+1 = since the measure in (|2.2|) is invariant under 
M -> —M. 

It can be easily seen that F(x) satisfies 

x 4 F 2 -(g + x 2 )F + g + x 2 + gx 2 U 2 = 0. (4.35) 

Hence, the generating functional has the form 



x 2 



+ g- \/(x 2 + g) 2 - Ax A (g + x 2 + gx 2 Tl 2 ) 
2x~ 4 



(4.36) 



Therefore, the Green functions n n are expressed in terms of n 2 , i.e. to 
know the Green series it is necessary to write down an equation for n 2 . This 
fact was discovered in |13j . One must mention that the approach based on 
the planar Sc winger-Dyson equations is not self-sufficient: it does not give 
such an equation on n 2 . One can write down the required equation within 
the approach proposed in [Ij. For the action ()4.33|) this equation looks like 

27g 2 U 2 2 + (1 + 18#)n 2 - 1 - 160 = 0. (4.37) 

2. Planar parquet generating functional 
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2r i-n2r ( 4 - 38 ) 




As it was said in section |21 the planar parquet approximation takes into 
account only a subset of all planar diagrams and is defined as an approxima- 
tive solution of the Schwinger-Dyson equations. In this limit the Schwinger- 
Dyson equations are reduced to a certain set of equations on the "basic" 
Green functions, i.e. Il 2 , T 3 , T 4 . For the action ()4.33|) this set looks like 

n 2 = 1 - 2gU 2 - g uir 4 
T 4 = -g + H + V 
H = -gU 2 T 4 + VI[ 2 2 T 4 
V = -gU 2 2 T 4 + HIl 2 T 4 

or 

= 1 - 2gU 2 - gU^Ti 

2^n|T4_ (4.39) 

^ + nir 4 -i 

The higher Green functions are constructed in terms of n 2 and T 4 . 

The first equation is nothing else but the Schwinger-Dyson equation 
(|4.34|) for n = 2. It can be easily seen since the full 4-point function II4 
in the planar parquet limit is 

n 4 = 2U 2 2 + n^r 4 , 

so 

n 2 = 1 - <?n 4 . 

Hence, the following equation on Il 2 derived from (|4.39p 

g 3 U 6 2 + g 2 U 5 2 + 5g 2 Tll + 5gU 3 2 + (1 - 5g)U 2 2 - 2U 2 + 1 = (4.40) 

is the required equation on U 2 . 
3. Parquet Higgs model 
Consider the following action 

S = -Tr M 2 + ^^Tr M 3 + — Tr M 4 . (4.41) 
2 3^ 4iV 1 ; 



The Schwinger-Dyson equations in this case are 



n-2 



n n + An n+1 + gU n+2 = ^n-i-2, n>2. (4.42) 



i=0 
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The functional and its equation are 



oo 



F(x) = ^H 
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(4.43) 



n=0 



x 4 F 2 - (x 2 + \x + g)F + x 2 + \x + g+ 



Tltx^x 2 + Xx + g) + Ii 2 x 2 {\x + g) + U 3 gx 3 = (4.44) 



Thus, the generating functional can be calculated in terms of 3 arbitrary 
constants ELl, n 2 , n 3 . 

The Higgs model in the vicinity of the true vacuum is the same as (I4.41j) 
with A = X(g). 

The planar parquet set of equations is given in Section It can be easily 
seen that the planar parquet approximation gives only two equations. To get 
the system closed one proposes the following equation on tadpole diagrams 



Hence, the planar parquet approximation gives the exact solution in terms 
of the generating functional (J4.43|) together with the set of equations on n l5 



In the action was considered in the planar limit by means of the 
generating functional technique. It was said that the arbitrary parameters 
(boundary conditions) can be derived by careful studying the holomorphic 
properties of the generating functional. 

5 Conclusion 

We have considered two different approaches to the matrix Higgs model. We 
have got the following results. 

1. The planar two-point functions for the symmetric one-cut f!3.23|) . non- 
symmetric one-cut cases (13.26)) and (|3.29j) have the following asymptotics as 

g -»• oo 



rii = An 3 . 



(4.45) 



n 2 , n 3 . 
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and in the two-cut case the planar two-point function (J3.32)) is 



2. The asymptotical expressions for the parquet two-point functions in 
the false vacuum (|2.7|) . (|2.9|) are 

= -l-12g + o(g), 
D® = -l-8g + o{g), 

4g 

and in the true vacuum the two-point function has the following asymptotic 

(EH 

7-^(4) 115 

2 % 2 4 v 7 
Hence, -D^, two-point Green function computed within the planar par- 

(3) 

quet approach in the true vacuum, coincides with G 2 , two-point Green func- 
tion computed within the planar approximation also in the true vacuum. 

As it was mentioned the parquet planar approach gives a very good agree- 
ment with the exact results in the one-cut case for the case of the positive 
mass square. In the case of the matrix Higgs model the situation seems more 
delicate because of the its multi-phase structure. Nevertheless, the planar 
parquet approximation leads to a rather good agreement at least in the small 
coupling limit. 

We have shown that in the planar parquet approximation it is possible 
to construct the generating functional for the Green functions. Besides, we 
have shown that the generating functional in the planar limit depends on 
the approximation chosen and, therefore, can be restricted to a subset of 
diagrams just by modifying initial conditions of the system. 
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